Entanglement generation through an open quantum dot: an exact approach 
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We analytically study entanglement generation through an open quantum dot system described 
by the two-lead Anderson model. We exactly obtain the transition rate between the non-entangled 
incident state in one lead and the outgoing spin-singlet state in the other lead. In the cotunneling 
process, only the spin-singlet state can transmit. To discuss such an entanglement property in the 
open quantum system, we construct the exact two-electron scattering state of the Anderson model. 
It is striking that the scattering state contains spin-singlet bound states induced by the Coulomb 
interaction. The bound state describes the scattering process in which the set of momenta is not 
conserved and hence it is not in the form of a Bethe eigenstate. 
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Introduction: We present a new exact approach to elec- 
tron entanglement generation in an open quantum sys- 
tem. Entanglement has attracted much attention in wide 
range of physics: it is a resource for quantum informa- 
tion processing and provides a new insight into quan- 
tum phase transitions in statistical physics [1, 2]. In 
most studies, the entanglement properties are discussed 
in closed systems in equilibrium. In order to study ma- 
nipulation of entanglement, however, we need to consider 
an open system out of equilibrium. Entanglement gener- 
ation using electrons in mesoscopic structures has been 
proposed recently [3-8]. In Refs. [3, 4, 8], in particular, 
devices are connected to reservoirs, electrons enter the 
device from the reservoirs, and interactions (the Coulomb 
interaction as well as the interaction between electron 
spin and nuclear spin) are essential for the entanglement 
generation. It is our purpose to discuss the entanglement 
generation in such a situation through an exact solution 
of scattering theory. 

In this Letter, we obtain an exact result for entangle- 
ment property of transported electrons of the two-lead 
Anderson model. The Anderson model is a fundamental 
model describing the electron transport through a quan- 
tum dot as illustrated in Fig. 1(a). It consists of one 
energy level (the quantum dot) in which two occupying 
electrons with opposite spins interact with each other 
(the Hubbard interaction), and two leads of noninteract- 
ing electrons each of which lead is coupled to the dot. We 
calculate the transition rate from the non-entangled inci- 
dent state with momenta k\ and hi on the lead 1 to the 
singlet and triplet states with momenta q\ and qi on the 
lead 2. In the scattering process which conserves the set 
of momenta as in Fig. 1(b), both the triplet and the sin- 
glet components of the incident state can be transmitted 
to the lead 2. On the other hand, in the cotunneling pro- 
cess which conserves the total energy but not the set of 
momenta as in Fig. 1(c), only the singlet component can 
be transmitted and the triplet component is filtered out. 
We clarify this mechanism by calculating the transition 
rates exactly, which is the main achievement of our ap- 
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FIG. 1: (a) A schematic diagram of the Anderson model, (b) 
A scattering process which conserves the set of momenta, (c) 
A cotunneling process. 



proach. The mechanism of the entanglement generation 
was first proposed in Ref . [3] ; the lowest order of our re- 
sult reproduces their perturbative result. An interaction- 
induced orbital entanglement property has been also dis- 
cussed in a quantum-dot system [9] . 

For the above purpose, we obtain the exact solution 
of the two-electron scattering state. A remarkable point 
of the state is that it contains a two-body singlet bound 
state. The bound state is induced in the cotunneling pro- 
cess (Fig. 1(c)) by the Hubbard interaction on the quan- 
tum dot. A many-body eigenstate of the closed Anderson 
model can be constructed by the Bethe Ansatz [10-13]. 
In contrast, our exact scattering state is a many-body 
eigenstate of the open Anderson model and essentially 
different from the Bethe eigenstate. A similar bound 
state is also discussed in Ref. [14] , where the two-electron 
scattering matrix has been constructed exactly in the An- 
derson model. While their study is focused on the asymp- 
totic states of electrons which lie far from the quantum 
dot, our exact solution describes electron states both in- 
side and around the quantum dot. 

Model and Result: The Hamiltonian of the Anderson 
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model is defined as H = Hq + Hi where 

/l d 

+ ^2 e d n da + Un d ^n dl , 



(1) 



H i = 7i E E ( c l(o)^ + 4*(o)) . (2) 



<r=T,i £=1,2 



Here (or) {ci„{x)) denotes the creation (annihilation) 
operator of an electron with spin <r(=T, [) on the lead I (= 
1,2). The dispersion relation is linearized to be e(k) — 
Vpk in each lead. Hereafter we set v F = 1. The operator 
d\ (d a ) represents the creation (annihilation) operator 
on the energy level of a quantum dot with energy e<j, and 
nda = d\d a . The parameter t represents the coupling 
between each lead and the dot. When the energy level is 
occupied by two electrons with opposite spins, they feel 
the Coulomb repulsion energy U > 0. 

We consider the situation studied in Ref. [3]. Let 
| k\ , k 2 ; 1) be the non-entangled incident state with mo- 
menta k\ and k 2 on the lead 1 defined by 



|fci,fc 2 ;l) =cl feiT 4- 2 j0) 



(3) 



where c\ k(T (cg kcr ) denotes the creation (annihilation) op- 
erator of an electron with momentum k and spin a on the 
lead £, and |0) denotes the vacuum state. We also define 
a triplet state \q\, q 2 ; 2, +) and a singlet state \q\, q 2 ; 2, — ) 
with momenta q\ and q 2 on the lead 2, 

l<?i>?2,2,±) = (4 qi Ai2l ±c 2 qi i4 q ^) |0>, (4) 

which are used as outgoing states. We calculate the tran- 
sition amplitude between these states: 



( qi ,q 2 ; 2,±\T{E k )\k u k 2 ; l)S(E k - E q ), 



(5) 



where E k — fci + k 2 , E q = q\ + q 2 and T(E) represents 
the transition matrix recursively defined by 



T{E)^H 1 +H 1 - 
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-T(E). 



E - Ho + i0 

We obtain new exact results for the transition rate: 

(gi, (ft; 2, +\T(E k )\k u k 2 ; \)5{E k - E q ) 
t 2 e kl e k2 



(6) 



4V2i 



(6(ki - qi)5(k 2 - q 2 ) - 5(k x - q 2 )S(k 2 - q^) 



(7) 



{q u q 2 ;2,-\T{E k )\k u k 2 ;l)5{E k -E q ) 



t 2 e kl e k2 
4V2i 



(<J(A;i - qi)S(k 2 - q 2 ) + Sih - q 2 )S(k 2 - 



where e k is defined in Eq. (19) below. Expanding them 
in the lowest order of t, we reproduce the perturbative 
result obtained in Ref. [3]. 

Equation (7) and the first term in Eq. (8) represent 
the contributions from the scattering process which con- 
serves the set of momenta {k\, k 2 } = {qi, q 2 } as shown 
in Fig. 1(b). Note that they vanish if ki and qj satisfy 
the condition of the cotunneling process (Fig. 1(c)), 



E k = E q , {fci, k 2 } ^ {qx, q 2 }. 



(9) 



In contrast, Eq. (8) has an additional term which re- 
mains non-zero for U > under the condition (9). The 
contribution appears only in the transition into the sin- 
glet state (Eq. (8)), not in the transition into the triplet 
state (Eq. (7)). In other words, we will observe only sin- 
glet states if we wait for outgoing electrons on the lead 2 
under the condition (9). 

Figure 2 shows the dependence of the transition rate 
2n\(qi, q 2 ; 2, — \T(E k )\k\, k 2 ; 1)| 2 on the interaction en- 
ergy U for t < E k (Fig. 2(a)) and t > E k (Fig. 2(b)). 
The solid lines indicate our exact result and the dashed 
lines represent the perturbative result in the lowest order 
of t obtained in Ref. [3]. The perturbative one, being 
divergent, fails when U ~ E k even for t < E k . 

Exact two- electron scattering state: In order to discuss 
the entanglement generation in such an open quantum 
dot system, what kind of quantum states should we focus 
on? What is the origin of the second term in Eq. (8)? 

We find that the two-electron scattering state \cj>) given 
by the solution of the Lippmann-Schwinger equation, 



\ki,k 2 ;l) + 



1 



-HM) (10) 



E k -H + 10 

plays an important role. In fact we have the relation 



T(E k )\k 1 ,k 2 ;l)=H 1 \ 



(11) 



which enables us to calculate the transition ampli- 
tude (5). 

In the remaining part of this paper, we construct the 
exact solution of \<f>). In Ref. [15], we obtain the exact 
solution of many-electron scattering scattering states in 
a spinless model. Here we apply the same technique to 
the two-lead Anderson model. The importance of \</)) 
results from the fact that the system is an open quantum 
system. It is suitable for the study of transport properties 
in the open quantum system, but the exact solution of 
such a state has not been discussed very much. As we 
will explain later, our exact solution has a remarkable 
property that the Bethe eigenstate does not have. 

By applying the transformation c ea (x) = 
(c la (x) + c 2rT (x))/V2, c oa (x) = (c la (x) - c 2l7 (x))/V2, 
we decompose the Hamiltonian into H = H e + H a , 
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FIG. 2: The transition rate 27r|(«i, q 2 ; 2, -\T(E k )\ki, k 2 ; 1)| 2 
obtained from our exact solution (solid line) and a per- 
turbative one (dashed line) for = 0, fci = k 2 — 0.1, 
gi = 0.05, q 2 = 0.15, E k = E q = 0.2 with (a) t = 0.1 (< E k ) 
and (b) t = 0.3 (> £?*)• 



where 



<T=T,1 



1 d 
rf.r 



c eC r(a;) + e<in dc 



£ t (4,(0)^ + 4(^(0)) 



(12) 
(13) 



Note that the odd part if is completely decoupled from 
the even part H e . Thus the two-lead Anderson model (1) 
and (2) can be transformed to the one-lead Anderson 
model (12) with the free part (13). Let \kx, k%; s) a p, 
{a, (3) — (e, e), (e, o), (o, o), s = +, — denote the scatter- 
ing eigenstate with the incident state being a two-electron 
plain- wave state 



l^^2;±)S = ^(4 lT 4 fai ±4 fai ct fc2T )|o). (14) 



€ 2 ;s> 



The scattering state \<j>) is expressed as [15] 
\<t>) = T/s(E (l^i' fc 2; s) ee + |/ci,fc 2 ; s) eo + |fci,/c 2 

* s=± 

- |fe 2 ,fci;+) eo + |A;2,fci;-) eo J, (15) 

because substituting the incident states (14) for all states 
on the right-hand side gives the incident state (3). 

Before giving the two-electron scattering eigenstate, we 
first mention the one-electron eigenstate \ka) e / defined 
as a solution of the Schrodinger equation H e / \ka) e i a — 
k\ko) e / . The one-electron eigenstate can be expressed 
as follows: 

Me = ^ dxg k (x)cl(x) + e k d£) |0), (16) 
\ka) = ( dxh k (x)cl a (x)\0) : (17) 



where the eigenfunctions are given by 
1 



9k(x) 



/2tt 



ek 



e lkx (0(-x)+6(x)^ 



1 



t 



h k (x) 



2irk-e d + it 2 /2 : 
1 



Akx 



/2tt 



(18) 
(19) 
(20) 



Here 9(x) is the step function and e k /e1 in Eq. (18) rep- 
resents the phase factor due to the scattering by the dot. 

Next we consider the two-electron scattering state, 
which is written in the form 

\ki,k 2 ]S = ±) a p 

dxe 3 a/3 (x;k l ,k 2 ) (c f ^(x)d\ ±ct /l (x)d\ S ) 
+ f^(k 1 ,k 2 )d\d^\0). (21) 

Here 7 = e for (a, [3) — (e, e), and 7 = o otherwise. From 
Eqs. (11) and (15), we find that the transition rate (5) is 
expressed in terms of the eigenfunctions: 

{q 1 ,q 2 ;2,±\T(E k )\k 1 ,k2]l}S(E k - E q ) 
t f°° 

= 5(E k -E q )—= dx(h qi (x) Th q2 (x)) 

8V27T J-00 

x (ef e {x;k 1 ,k 2 ) - ef (x;k u k 2 ) ± ef {x;k 2 ,k 1 )). (22) 

We construct the eigenfunctions specifically e s a Ax) for 
all cases of (a, (3) to calculate the transition rate (22). 
First of all, let us consider the case a = (3 = e, s = 
— . Only in this case, the eigenfunctions depend on the 
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Coulomb interaction U. The functions with total energy 
Ek = ki + k 2 are obtained by solving the two-electron 
Schodinger equation 



(<9i + d 2 ) - E k g ee (xi,x 2 ) 



+ t (S{xi)e ee {x 2 ) + 5(x 2 )e ee (x 1 )) = 0, 



(23) 



~ + td - E^j e-(x) + tg-{x, 0) + tS(x)f~ = 0, 

(24) 

(2e d + U- E k )f~ + 2te" (0) = (25) 
under the condition 

g-(x u x 2 ) = —L- f e ^ klXl+k2X2) + e i( ^ klX2+k2Xl A (26) 
2y27r ^ ' 

for x\,x 2 < 0. The result is as follows: 

g- e (x x ,x 2 ) = (^Sfci ( X Q i)9k 2 {xq 2 ) 



+ 



u f^ZX b -^ " a ' z(x ^ )eM )- (27) 



e ee ( x ) = (ek 1 9k 2 (x) +e k2 gk 1 (x)) 

V2itUe kl e k2 j EkXZ{xl 



+ 



fee = V2e kl e k2 - 



2e d + U-E k - if 2 

V2Ue kl e k 



2e d + U-E- it 2 



(28) 
(29) 



where Q = (Qi,Q 2 ) is a permutation of (1,2), iy 



Xi — xj and 



Z{x) = e^ d+t2 ^ x d(-x). 



(30) 



It is obvious that Eq. (27) satisfies the condition (26). 
Note that the condition (26) corresponds to the plane- 
wave incident state (14). We confirm that our solu- 
tion (27)-(29) of the Schrodinger equation (23)-(25) sat- 
isfies the Lippmann-Schwinger equation, 

\ki,k 2 ;-) ee = \h,k 2 ;-)^ + — 1 - ffi|fci, k 2 ; -) ee . 

tjk — no + lU 

(31) 

Remarkable is that Eqs. (27) and (28) contain the two- 
body bound state Z(x) in Eq. (30). The range of binding 
is t~ 2 , which itself is independent of U. Note the follow- 
ing properties: (i) It describes a scattering process that 
does not conserve the momentum set as the cotunnel- 
ing process (9). In contrast, the first term of Eq. (27) 
represents the direct (fci = q\,k 2 = q 2 ) and exchange 
(fci = q 2 ,k 2 = qi) processes, (ii) It is induced by the 
Coulomb interaction U; It vanishes for U — 0. (iii) It 
only appears in the case a — (3 = e and s = — since f^g 
in Eq. (21) vanishes in the other cases and thereby the 



eigenfunctions do not depend on U. From these proper- 
ties, we find that it is this bound state that produces the 
second term of Eq. (8). 

This solution is essentially different from the Bcthc 
eigenstate [10-13]. In the Bcthc ansatz, we would sup- 
pose that g~ e (xi,x 2 ) is in the following form: 

g~ e (x u x 2 ) = ^A Q g kl {x Ql )g k2 {xQ 2 ). (32) 
Q 

Here Aq depends on ki and the parameters of the model 
(td, t and U). Note that the solution is characterized by 
the fixed set {k\,k 2 }. Although fc» could be a complex 
number, it is obvious that Eq. (32) cannot take the form 
of Eq. (27). 

In the different cases of a, (3, s from a = (3 = e, s = — , 
the scattering states are interaction-free eigenstates. In 
particular, e s aB {x) are obtained as follows: 



ete(x) = (-e kl g k2 (x) + e k2 g kl (x)) . 



(33) 



1 



efoix) = T^e kl h k2 (x), ef (x) = 0. (34) 

The scattering state \<j>) constructed from Eq. (15) with 
these eigenfunctions is shown to satisfy the Lippmann- 
Schwinger equation (10). Using all these exact solutions 
in Eq. (22), we finally arrive at the desired results (7) 
and (8). 

Conclusion: In this Letter we have constructed the 
exact two-electron scattering state and discussed its en- 
tanglement property. For the exact calculation of the 
quantities (7) and (8), our solution is essential. We have 
clarified that the electron transport through the quantum 
dot has a potential advantage for entanglement genera- 
tion. The construction of many-electron scattering states 
in other more complicated models and the understand- 
ing of their entanglement property are interesting issues 
in the future. 

Our exact scattering state could also be a powerful 
tool to understand non-equilibrium electron transport 
with finite bias voltage. In a spinless model, the ana- 
lytic approach for this topic has been proposed [16-19] 
and the importance of many-electron scattering states 
on the noncquilibrium current has been pointed out re- 
cently [15, 20, 21]. In particular, our approach in Rcf. [15] 
succeeded in obtaining nonperturbative result on nonlin- 
ear current-voltage characteristics. This suggests that 
our exact many-electron scattering state yields impor- 
tant information about the Kondo effect out of equilib- 
rium [22-25]. 
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